Thermal rectification in nonlinear quantum circuits 
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We present a theoretical study of radiative heat transport in nonlinear solid-state quantum cir- 
cuits. We give a detailed account of heat rectification effects, i.e. the asymmetry of heat current 
with respect to a reversal of the thermal gradient, in a system consisting of two reservoirs at finite 
temperatures coupled through a nonlinear resonator. We suggest an experimentally feasible super- 
conducting circuit employing the Josephson nonlinearity to realize a controllable low temperature 
heat rectifier with a maximal asymmetry of the order of 10%. We also discover a parameter regime 
where the rectification changes sign as a function of temperature. 



I. INTRODUCTION 

Heat transport in nanoscale structures has become an 
active and rapidly growing research area. Progress in 
experimental methods has enabled the study of funda- 
mental issues, and lately the field has seen major break- 
throughs, such as the measurement of quantized heat 
transport and manipulation of thermal currents using 
external control fields.—^ In solid-state systems electron- 
electron and electron-phonon scattering are the most im- 
portant channels for small systems to exchange energy 
with the environment. However, recently it was under- 
stood that at low temperatures one needs to take into 
account the radiative channel which becomes the domi- 
nant relaxation method in mesoscopic samples below the 
phonon-photon crossover i^i^ 

In this paper we study rectification effects in thermal 
transport mediated by electromagnetic fluctuations in 
solid-state nanostructures. In a two-terminal geometry 
a finite rectification means that heat current is not sim- 
ply reversed when the thermal gradient changes sign, but 
also the absolute magnitude of the current changes. We 
define the rectification TZ as 



TZ = (J+ — J-)/ max{J+, J_}, 



(1) 



where J+ and J_ are the magnitudes of the heat cur- 
rents in forward and reverse bias configurations, respec- 
tively (see Fig. [T]). Previously rectification has been 
shown to take place in systems where a classical^i^ or 
quantized^ii^ nonlinear chain is coupled asymmetrically 
to linear reservoirs, when nonlinear reservoirs are cou- 
pled through a harmonic oscillator, or in hybrid quan- 
tum junctions. Here we demonstrate rectification in 
a fully quantum-mechanical and experimentally realiz- 
able model where photon-mediated heat current flows 
between two linear reservoirs coupled asymmetrically to 
a nonlinear resonator. 

Our analysis is based on a nonequilibrium Green's 
function method developed in Ref. [isl . and the nonlinear 
transport problem is solved with a self-consistent Hartree 
approximation. Rectification is studied as a function of 



2high [ 



forward bias 




Ml Mr 



2J-0W "Xf 



reverse bias 
J_ 



^high 




Ml M, 



FIG. 1: In a heat transport experiment, thermal energy flows 
from a hot reservoir (temperature Thigh) to a cold reservoir 
(temperature Tiow). To obtain the rectification TZ one must 
measure the current for both thermal bias directions. In our 
model the heat is transported by inductive magnetic coupling 
between the reservoirs and the central nonlinear resonator. 



the operating temperatures, reservoir coupling strengths 
and admittances, and the strength of the nonlinearity. 
We also propose a concrete setup based on a super- 
conducting quantum interference device (SQUID) where 
the rectification effects can be realized with current ex- 
perimental technology at sub-Kelvin temperatures. A 
similar circuit, operated in the linear regime, was em- 
ployed in the pioneering experiment demonstrating pho- 
tonic heat transport;^ By adjusting the external magnetic 
flux through the circuit it is possible to tune the rec- 
tification continuously between zero and the maximum 
value. Using realistic parameters we find a rectification 
of over 10%, and identify a regime where TZ changes sign 
as a function of temperature. Experimentally rectifica- 
tion has been observed in phonon transport through a 
nanotubei^ at room temperature with TZ = 7% and in 
electron transport through a quantum dot^^ at 80 mK 
with TZ up to 10%. 



II. MODEL 

The thermal transport setup is depicted in Fig. [TJ It 
consists of two linear reservoir circuits with admittances 
Yl{uj) and Yii{iv). Temperatures of the left and right 
reservoirs are Thigh and Tiow < Ihigh in the forward bias 
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setting, and vice versa for reverse bias. We assume that 
heat can flow between the reservoirs only through a medi- 
ating nonlinear resonator circuit. The couplings between 
the reservoirs and the resonator are taken to be induc- 
tive with mutual inductances and M^. Using the 
Caldeira-Leggett mapping between linear admittances 
and bosonic reservoir modes the total Hamiltonian takes 
the form H — Hl + Hr + Hm + He , where the middle 
circuit and reservoir terms are 



1, 



he 



Hl/r — 



E 

jeL/R 



and the inductive coupling term is 



Hc^li MLiL + Mrir 



(2) 



(3) 



(4) 



which involves the current operators for the central de- 
vice / and for the reservoirs il/r — J2j£L/R dji^J + '^J)' 
respectively. The electric current operator for the cen- 
tral device can be expressed as / = Io{b + b'') with 



Iq = ^Jhui^jlL and ujq = l/V LC where L and C are 
the linear inductance and capacitance of the resonator; 
6, w and reservoir operators are bosonic creation and an- 
nihilation operators, [6,6^] = 1. The nonlinearity of the 
central circuit is characterized by the the second term 
in Eq. ([2]), corresponding to a quartic potential whose 
strength is controlled by the parameter e. It must be 
emphasized that Eq. (jU has a generic bilinear form and 
therefore our results are relevant for other types of sys- 
tems beyond the studied realization. 

The basis of our analysis is provided by the Meir- 
Wingreen formula for the heat curren t ^^i"'^^ 

J = r {2 [S,{lo) - Si{-u:)] MYL{oj)]nL{uj) 

Jo 



2Si{-uj)Re[YLiiu)]}. 



(5) 



Here nL(uj) is the Bose function of the left reservoir and 
Si{uj) = J^^dte"^^*~*'\i{t)i{t')) is the current noise 
power of the central circuit. The admittances Yi^/ji{uj) 
are related to the current correlation functions of the 
free reservoirs. In the absence of the nonlinear term 
(e — 0) the transport problem can be solved exactly for 
arbitrary couplings and reservoir admittances,''^- No rec- 
tification takes place in this regime. In the following we 
solve the nonlinear transport problem in a self-consistent 
Hartree approximation, which is expected to be accu- 
rate for small values of the nonlinearity. This approach 
does not fully account for the correlation effects due to 
the interplay of nonlinearity and tunneling which are po- 
tentially important in the ultra-low temperature regime 
Thigh, 2]ow ^ ^o/^s- However, analogously to interact- 
ing electron transport problems, the mean-field approach 



is accurate in the sequential tunneling regime when the 
temperatures are of the order of hw^/kB- 

As a first step we approximate the resonator Hamilto- 
nian as 



Hm 



huJoib'^b+h + ihe^ib'^ + bf 



(6) 



where we have used (6+ + h)"^ « 6>I>(6t + bf with the 
mean field $ ((6^ + b)"^)- Here the factor 6 is the 
number ways two operators can be picked from a set of 
four. By performing a diagrammatic expansion of the res- 
onator Green's function one can show that this procedure 
is identical to the self-consistent Hartree approximation. 
Because Eq. (|6]) is now quadratic in bosonic operators, 
it is possible to bring it to diagonal form by a canonical 
transformation. However, now we have the added com- 
plication of an a priori unknown mean field, which has to 
be evaluated self-consistently in a nonequilibrium state. 
The transformed Hamiltonian and current operators are 



H 



M 



h{b + b^), (7) 



where a)o = t^o^ 1 + ^ and /q = y^^/o- Thus the ef- 
fect of the nonlinear term is incorporated by a mean-field 
dependent renormalization of the resonance frequency 
of the oscillator and its current operator. For further 
development it is convenient to introduce the correla- 
tion functions {i{t)i{t')Y = -ie{t - t'){[i{t),i{t')]) and 
{i{t)i{t'))< = -i{i{t')i{t)). A nonequilibrium equation- 
of-motion analysis^^, similar to the one presented in 
Ref. [l^ reveals that the current correlators are given 
by 



(//)<(c.) 



V2|(/7yMpi]<H, 



(8) 
(9) 



where (//)q(cli) = 2l'^ujo/ {oj'^—Cjq) is the retarded Green's 
function of the uncoupled oscillator. The self-energies 



S'-(C.) : 

E<H = - 



MIYl{u) + MIYr{u;) , (10) 



MlRc[YL{Lo)]nL{Lo) 



h 

+ MlRe[YR{Lo)]nR{uj) 



(11) 



take into account the presence of reservoirs. Further- 
more, the mean field $ is related to the lesser correlator 
via 

$ = ((6t + if) = -l-^ 1^ (12) 



Equations (l8|)- (|T2|) form a closed set of equations which 
needs to be solved to find the current correlation func- 
tions. The self-consistent solution proceeds by making 
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an initial guess for the mean field, calculating the cor- 
relation function ^ corresponding to the initial value 
and calculating the updated value of the mean field 
by evaluating the integral in Eq. (jl2p . The procedure 
is repeated until convergence is achieved. The current 
noise then follows immediately from the lesser function 
Si{uj) = — Im( //)< (— (jj) which yields the heat current 
after evaluating Eq. ^ . In the case of a vanishing non- 
linearity (e — 0) this procedure recovers the exact solu- 
tion of the linear problem. To facilitate the analysis of 
the rectifying mechanism we note that with the help of 
Eqs. ©-(H^) we can write Eq. ^ in the form 

r dco Ahu^MlM]iRe[YL]Re[Yi^]inL - ur) ^^^^ 
Jo 27r + iuj[MlYL + Aflli?] 1^ 

where F{uj) = h{uj'^ - ~ I2ujf)e^) / {2I^loq) . The fre- 
quency dependence of YL/Ri''^) and nL/ulu!) has been 
suppressed for brevity. 

For numerical calculations explicit expressions for the 
admittances Y^/.rC'^) are needed. Here we assume that 
the reservoir circuits effectively consist of a resistor, a ca- 
pacitor, and an inductor in series, resulting in ^^/^(aj) = 

- 'QL/nizjf;^ - ^)]-\ where R^/n, Ql/r, 
and iui^/n are the resistance, quality factor and resonance 
frequency of the left and right reservoir, respectively. The 
behavior of the system is now uniquely determined by 
nine dimensionless parameters: e/wo, fcsTiow/high/^Oj 
^L/Rr^o /^^L/Ri Ql/Ri and uji^fji/ujQ. Rectification can 
then be calculated from Eq. ([T]), by computing the for- 
ward and reverse bias currents, with the above 
prescription. 



III. RESULTS 

Let us illustrate some generic features of the model 
with the simple setup of two purely dissipative reser- 
voirs, Ql = Qr = 0, in which case the frequencies ujl 
and lor are irrelevant. In Fig. [5] we plot the rectification 
against three different variables. First, from Fig. [H^a) 
we see that already at quite small values of nonlinearity, 
e ~ 0.1 Wqj the rectification has essentially reached its 
maximum. Such values for e are well within the regime 
of validity of our approximations and should also be eas- 
ily achieved in the experimental setup proposed below. 
Next, Fig. EJb) exemplifies a very generic feature: hav- 
ing M1/Rl < M'^/Rr tends to produce J+ > and 
vice versa. Finally Fig. ^c) shows that the rectifica- 
tion increases logarithmically with the temperature ratio 
Thigh/Tlow Therefore, to see an appreciable effect, the 
temperature difference Thigh — Tiow should be of the same 
order of magnitude as the temperatures themselves. 

For purely resistive reservoirs maximal value for the 
rectification is about 2% (Fig.[21[a)). Larger values can be 
obtained by adding a reactive part to one of the reservoir 
circuits. Then, as Fig. [3] shows, TZ can be made an order 




FIG. 2: Rectification with purely resistive reservoirs, Ql = 
Qr = 0, as a function of (a) nonlinearity e, (b) coupling 
Ml, and (c) temperature ratio Thigh /Tiow. In all panels we 
have e/cJo = 0.07, fcsThigh/fiti^o = 0.2, fcsTiow/ftwo = 0.1, 
MIIq/HRl = 0.2, and M^Io/hRR = 1, except for the vari- 
able on the horizontal axis. In panel (c) Thigh is varied. 

of magnitude higher. The inset shows the current J+, 
normalized with respect to the universal single-channel 

maximum heat current Jmax = ^(T^igh - ^low)-'^ Ac- 
cording to Fig. [31 the highest values for TZ are obtained 
for high temperatures, where J+ tends to zero. High rec- 
tification and large current are thus competing effects, 
and the optimal operating point depends on the experi- 
mental constraints. In any case, it is possible to obtain a 
rectification of ~ 5% with J ~ 0.1 Jmax and up to ~ 15% 

with J ^ 0.01 Tniax- 




FIG. 3: Rectification with one reactive reservoir {Ql = 0.1). 
Here Thigh /Tow = 2 and the different curves correspond to 
lol/ujo = 0.2 (solid), 0.1 (dashed), 0.05 (dash-dotted), 0.02 
(dotted). Other parameters as in Fig. [5] 

From Fig. [3] we also see that decreasing ujl in- 
creases rectification, so both small lul and the condition 
M1/Rl < M'^/Rr favor the direction J+ > J_. We can 
also combine these two trends in an opposing manner by 
making uj^ large. This way one can produce a system 
where the direction of rectification changes as a function 
of temperature. From Fig.|4]we see that in a system with 
a high-frequency reservoir (here ujl = 10 wq), TZ is posi- 
tive when both temperatures are below Huo/kB, but at 
higher temperatures the same device produces a negative 
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FIG. 4: Change of rectification sign as a function of operat- 
ing temperature. Here the lower temperature is kept fixed, 
with values ksTio^ /hujQ = 0.1 (solid), 0.3 (dashed), 1 (dash- 
dotted), 3 (dotted), and the higher temperature is varied. The 
left reservoir is reactive with Ql = 0.1 and ujl = lOtJo, and 
other parameters are as in Fig. [21 



TZ. In contrast to previous reportsSii^ on the rectification 
sign reversal, in our system only the reservoir tempera- 
tures need to be changed, not the device parameters. 

The direction of rectification can be understood as fol- 
lows. Equations ^ and show that the mean field 
$ is proportional to the self-energy E^. Due to the self- 
consistency loop the relationship between $ and S< is 
not actually linear, but in practice making larger will 
also increase $. Next, Eq. pT|) shows that S< is es- 
sentially the product of Bosc function n^/^(w) and the 
effective coupling strength a;A/|^^Re[YL/i{(a;)], summed 
over the two reservoirs. Because of this form, it follows 
that when comparing the forward and reverse bias set- 
tings, larger is obtained in the case when the more 
strongly coupled reservoir is hotter. As a consequence, 
the mean field $ is also larger when the more strongly 
coupled reservoir is hotter. 

To interpret physically the ^-dependence of the cur- 
rent J, we analyze separately the numerator and denom- 
inator of the integrand in Eq. p3p . The numerator is the 
product of the energy hw, Bose window 71^(0;) — nfi{iLj), 
and effective reservoir coupling strengths, as defined 
above. Thus it can be seen as a measure of the energy 
available for transport at the reservoirs. On the other 
hand, the denominator is due to the Green's function 
\{Iiy {u!)\'^ , giving the transmittance of the central cir- 
cuit. Equation shows that an increasing $ effectively 
increases the central circuit resonance frequency, thereby 
shifting the resonator transmission window to higher en- 
ergies. Because of the Bose functions, in most situations 
the numerator is smaller at higher energies and the total 
current decreases. But this is not always the case. It 
turns out that if at least one of the reservoirs is reactive 
with high resonance frequency (> luq) and the reservoir 
temperatures are high (> hujo/kB), the peak of the nu- 
merator is shifted to high enough energies so that an 



increasing $ produces on increasing J. In summary, ex- 
cept for the case of high-temperature and high-frequency 
reservoirs, larger current is obtained in the configuration 
where the more weakly coupled reservoir is hotter. This 
explains the sign of TZ in all our results. 



IV. EXPERIMENTAL REALIZATION 

For low operating temperatures, with Thigh, Tlow ap- 
proximately in the range 100 mK-1 K, the studied model 
can be realized by the setup shown in Fig.[5l The system 
consists of a superconducting loop containing a Joseph- 
son junction characterized by its Josephson energy Ej 
and shunt capacitance C. The loop itself is assumed to 
have a finite inductance dominating the potential land- 
scape. The Hamiltonian of the system is^^ 

HM^Ecf + EL{$-(j),f -Ej cos $, (14) 

where the charging and inductive energies are Eq = 
e^/2C, El — (/i/2e)^/2L, and (l)x denotes the exter- 
nal magnetic flux through the loop (in units of h/2e). 
The superconducting phase across the junction cj) and the 
charge at the capacitor q (in units of electron charge) are 
treated as conjugate observables [(j), q] = 2i. The charg- 
ing term can be thought of as the kinetic energy and 
the 0-dependent terms as an effective potential energy 
of a fictitious particle. In the following we assume that 
(j>x Ri TT and Ej < 2El so that the potential has a single 
minimum a,t (f) = 4>o, with (j)Q « tt. With these assump- 
tions the phase is bound close to the minimum so that we 
can approximate the potential accurately by expanding 
the cosine term to the 4th order: 

Hm Ecf + (El + ^Ej cos M^'^ - ^Ej cos cjjo 4>'^ 
= Ecf + E2$^ + Ei$\ (15) 

the second line defining the quantities E2 and E4. In 
general there should also be a (f)^ term, but with 0o ~ tt 
this is small. Further, within the mean-field approxi- 
mation one has cjr' ~ (p{(/)^), producing just a shift in 
the origin. Writing the charge and phase in terms of 
bosonic creation and annihilation operators we recover 
exactly Eq. ^ with parameters Hcuq = A\/EcE2 and 
he = The current operator of the circuit is given 

by / = /o(6 + St), where /q = AeiEcElY''^ /h. Thus, 
in the parameter regime E2 S> -B4 we have effectively 
realized our weakly nonlinear resonator model. 

As the above considerations show, varying the exter- 
nally applied field (j)x about tt moves the potential min- 
imum (pQ which in turn changes the values of the pa- 
rameters L0{) and e. In particular, e is maximized at 
4'x — fpa — and vanishes when — > tt ± 7r/2. Fig- 
ure [5] demonstrates the resulting continuous tuning of 
rectification performance. 
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FIG. 5: Rectification and heat current through the SQUID 
system as a function of the control flux (px- Parameter values 
at 0a: = TT as in Fig. [2] except Thigh ~ 2Tiow ~ OAhujo/kB- 
Inset: Schematic of the SQUID setup with two linear reser- 
voir circuits inductively coupled to a superconducting loop 
containing one Josephson junction. 



V. CONCLUSIONS 



In conclusion, we have analyzed heat rectification ef- 
fects in radiative heat transport through a nonlinear 
quantum resonator. This system is particularly inter- 
esting because it can be realized by an experimentally 
feasible superconducting circuit. The proposed system is 
operated in a low-temperature regime and can be con- 
trolled by by applying external magnetic fields. Despite 
its simplicity, the system is capable of producing a recti- 
fication of over 10%. We have shown that in a suitable 
parameter regime the direction of rectification changes as 
a function of temperature and given a physical explana- 
tion for the phenomena. 
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